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Abstract. We examine the spatial clustering of extremely red objects (EROs) found in a relatively large survey 
of 700 arcmin 2 , containing 400 galaxies with R — Ks > 5 to Ks = 19.2. A comoving correlation length ro — 12 ±3 
/i _1 Mpc is derived, under the assumption that the selection function is described by a passively evolving early- 
type galaxy population, with an effective redshift of z ~ 1.2. This correlation length is very similar to that of 
local L* elliptical galaxies implying, at face value, no significant clustering evolution in comoving coordinates 
of early-type galaxies to the limiting depth of our sample, z ~ 1.5. A rapidly evolving clustering bias can be 
designed to reproduce a null result; however, our data do not show the corresponding strong reduction in the 
average population density expected for consistency with underlying growth of the mass-function. We discuss our 
data in the context of recent ideas regarding bias evolution. 

The uncertainty we quote on ro accounts for the spikey redshift distribution expected along relatively narrow 
sightlines, which we quantify with detailed simulations. This is an improvement over the standard use of Limber's 
equation which, because of its implicit assumption of a smooth selection function, underestimates the true noise 
by a factor of « 3 for the parameters of our survey. We propose a general recipe for the analysis of angular 
clustering, suggesting that any measurement of the angular clustering amplitude, A, has an intrinsic additional 
uncertainty of a a/ A = VAC, where AC is the appropriate integral constraint. 
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cD; Galaxies: formation; Galaxies: fundamental parameters 
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former class: Broadhurst & Bouwens (2000), Moriondo, 
Cimatti & Daddi (2000) and Stiavelli & Treu (2000) have 
concluded from independent datasets that most (~ 70%) 
such objects have De Vaucouleurs profiles, with only about 
15% of them displaying irregular or disturbed morpholo- 
gies, expected for dusty starburst systems (the remain- 
ing 15% has a disk-like exponential profile). The existing 
spectroscopic results for single objects or for small sam- 
ples of EROs support this conclusion (Spinrad et al. 1997; 
Soifer et al. 1999; Cimatti et al. 1999, Liu et al. 2000). 
Spectroscopy of Mux-limited samples of K selected galax- 
ies generally suffer from incompleteness of the optically 
reddest galaxies but the weight of evidence is that most of 
the reddest objects have spectra consistent with early-type 
galaxies, up to the effective spectroscopic limit of z ~ 1.3 
for absorption-line work (Cohen et al. 1999; Eisenhardt 
et al. 1998; Cimatti 2001). These surveys also broadly 
show that the EROs of known redshift are in the range 



0.8, 



1.5, consistent with expectations for passively 



evolving ellipticals based on an extrapolation of the local 
luminosity function with passive evolution (Sect. |3j). 

Even if a few EROs have been identified as counter- 
parts of SCUBA sources (Smail et al. 1999; Gear et al. 
2000), SCUBA observations of complete ERO samples 
show no frequent detections (Mohan et al. 2001, in prepa- 
ration) reinforcing the idea that dusty HRlO-like objects 
(Cimatti et al. 1998, Dey et al. 1999, Andreani et al. 2000) 
are rare among EROs. It is therefore reasonable to con- 
clude that with K ,$ 19 EROs we are observing the clus- 
tering signal of predominantly distant early-type galaxies. 

Here we take the angular correlation function mea- 
surements of the EROs and a plausible estimate of their 
selection function to derive their spatial clustering am- 
plitude. The 3D correlation length of EROs should thus 
produce constraints on the evolution of the clustering of 
early-type galaxies, a population of objects which is likely 
to have formed in the highest amplitude perturbations and 
to be positively biased with respect to the general galaxy 
population and to the overall distribution of mass. The 
evolution of the correlation amplitude of early-type galax- 
ies is an observable which is independent of measurements 
of the evolution of their comoving number density, provid- 
ing therefore a complementary mean to examine the rate 
of evolution. Constraints on the density evolution of EROs 
have been sought previously in response to predictions of 
CDM models (Baugh et al. 1996, Kauffmann 1996). Early 
work based on small fields claimed observational evidence 
for a sharp decline in space density of early type galaxies, 
while more recent estimates based on deeper and more 
complete samples are consistent with a constant comov- 
ing density of this population up to at least z ~ 1.3, and 
imply a typical formation redshift of the stars in these 
galaxies not less than Zf ~ 2.5 assuming passive evolution 
(see Daddi, Cimatti & Renzini 2000 for a complete dis- 
cussion) . Hence it is very important to obtain information 
regarding the clustering evolution to independently ad- 
dress this important question. Here we analyze together 



both these questions of clustering and density evolution of 
EROs, with the largest complete sample of relevant data. 

In Sect, |2j we describe the standard Limber's equa- 
tion formalism. Sect. |3J derives order-of-magnitude results 
on the correlation length rg that confirm the interpreta- 
tion of EROs as high-z early type galaxies and justify 
our assumed redshift distribution. Sect. |j present the ro 
estimates based on Limber's equation. In Sect. |5j we dis- 
cuss the limitation of the standard approach and we use 
numerical simulations to find the definitive constraint on 
the correlation length, that differs significantly from the 
Limber's equation results. We discuss the general impli- 
cations of our findings for the clustering analysis on small 
areas. We than compare in Sect. Oour estimates of the cor- 
relation length of distant ellipticals to the measurements 
in the local universe and to theoretical models predictions. 
Our conclusions are presented in Sect. ^. 

All the scales quoted in the paper are given in comov- 
ing units. Three cosmological models have been consid- 
ered: a A-fiat universe (fl m = 0.3, tt\ = 0.7, h = 0.7), an 
open universe (f2 m = 0.3, fl\ — 0, h = 0.7) and an Sl-flat 
universe (0 OT = 1, il\ = 0, h = 0.5). Ho = lOOh km s _1 
Mpc -1 . 

2. The Limber's equation and the spatial 
correlation length 

The angular two point correlation function w(9) is re- 
lated to its real space analogous £(r) by Limber's equa- 
tion (Peebles 1980). In the case of small angles (9 <C 1), 
if both w and £ have power law shapes, writing £(z) = 
(r/ro(z)) _7 (with 7*0(2) being the comoving correlation 
length at redshift z, and r the comoving distance), the 
Limber's equation becomes: 

,„(ffS r- n{l-l)mS9{z)(dNldzfr,{zydz 

where dN/dz is the redshift selection function of the sam- 
ple, which in the limit of a large number of objects coin- 
cides with the observed redshift distribution. The function 
g(z) depends only on the cosmology: 

g{z) = (dx/dz^x^^Fix) (2) 

where x and F(x) arc defined with the metric: 

ds 2 = c 2 dt 2 - a 2 [dx 2 /F(x) 2 + x 2 {d9 2 + sin 2 9d<j) 2 )]. 

If we define: 



'O.eff 



then by using eq. [D and with w(9) = A9 1 7 we have: 



B=^ 



g(z)(dN/dz) 2 r a (zydz / / g(z)(dN/dz) 2 (3) 

(4) 
(5) 



T ((j-l)/2)Jg(z)(dN/dz) 2 dz 
T( 7 /2) [J(dN/dz)dz] 2 
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Thus, knowledge of a measured or assumed red- 
shift distribution allows us to relate the angular cluster- 
ing amplitudes A to the 3D correlation length r . o ff. If 
d 2 r (z) 7 jdz 2 is negligible in the relevant redshift range, 
then ro, e ff = r (z e ff) with: 



z e & 



zW(z)dz 



W(z) = g(z)(dN/dz) 2 / / g(z){dN/dz) 2 dz 



(6) 



(7) 



This inversion process provides a weighted estimate 
of ro over the probed redshift range. In the following, we 
will refer generically to ro as ro. c ff = ro(z e ff), bearing in 
mind, anyway, the effect of eqs. || and g. For consistency 
with D2000, where the clustering amplitudes were fitted 
with w(6) ex 9~°- s , we assume 7 = 1.8 for our analy- 
sis, consistent with most previous observations (see e.g. 
Roche & Eales 1999). In the following we will quote the 
values for the amplitude of the angular two-point correla- 
tion function as measured (or extrapolated) to 1 degree, 
i.e. A = A{1°). 

3. The effects of the selection function 

To derive 3D information from the angular correlation 
measurements a selection function for EROs must be sup- 
plied. As discussed in the introduction, strong evidence 
exist that the bulk of EROs is made of early type galax- 
ies. 

The strong angular clustering of EROs reported by 
D2000 and McCarthy et al. (2000), independently sup- 
ports this same conclusion. In fact, early- type galaxies 
in the local universe are known to prefer high density 
environments (Dressier 1980) and to have much larger 
correlation lengths than late-type galaxies (ro ^ 7-8 ver- 
sus r £ 5; Davis & Geller 1976; Giovanelli, Heynes & 
Chincarini 1986; Loveday et al. 1995) and than dusty 
starburst (see e.g. Saunders et al. 1992 for IRAS galax- 
ies that have r ~ 3.8 hr 1 Mpc). At higher (z ~ 1) red- 
shift even lower correlation lengths have been observed for 
star forming galaxies: Adelberger et al. (2000) find that 
balmer break (z ~ 1 star forming) galaxies have ro ^S 3 
hr 1 Mpc, while the blue starburst selected field galaxies 
at 0.8 < z < 1.5 have 1 < r < 2.5 h^ 1 Mpc (Carlberg et 
al. 1997, Hogget al. 2000). 

In Fig. [l] we calculate the width of simple top-hat red- 
shift distribution that can reproduce the ERO observed 
angular amplitude, as a function of ro. Requiring that 
EROs have a correlation length of ro ^, 3 hr 1 Mpc, typ- 
ically measured for z ~ 1 star forming galaxies, would 
imply a very narrow ERO redshift distribution of width 
Az rs 0.05 to reproduce the observed amplitude, which 
seems implausible. For any reasonably broad redshift dis- 
tribution a much larger ro is inferred, favoring the larger 
ro known for local early- type galaxies, but irreconcilable 
with the small correlation lengths observed for star form- 
ing galaxies at both low and high-z. 
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Fig. 1. The correlation length ro that reproduces the typical 
ERO clustering of A = 0.02 (D2000), as a function of the width 
of a top-hat redshift distribution centered at z — 1 (solid line) 
and z — 1.5 (dotted line), for the A-fiat cosmology. 



In the remainder of the paper the selection functions 
expected for distant early type galaxies will thus be used. 

3.1. Modeling the selection function 

We adopt models accounting only for passive evolution, 
appropriate for early-type galaxies, to estimate the ERO 
selection function. This is well justified for the present 
analysis that deals with the clustering of galaxies redder 
than R — Ks > 5, corresponding to z £ 0.8-0.9 for a pas- 
sively evolving L* elliptical, as several studies have shown 
that at least up to z ~ 1 the photometric and density 
evolution of the elliptical galaxies is consistent with pas- 
sive evolution with no number density evolution (Totani 
& Yoshii 1997; Franceschini et al. 1998; Schade et al. 1999; 
Im et al. 1999; Broadhurst & Bouwens 2000; Scodeggio & 
Silva 2000; Phillipps et al. 2000; Daddi, Cimatti & Renzini 
2000). 

For the passive evolution models adopted here, ellip- 
ticals form with a rapid burst (jsfr = 0.1 Gyr). The 
Salpeter IMF is assumed, with no dust reddening, and 
Z = Zq. The Bruzual & Chariot spectral synthesis mod- 
els (1993) in the 1997 version were used, with the Marzke 
et al. (1994) pure-ellipticals luminosity function for the 
normalization at z = 0. Daddi, Cimatti & Renzini (2000) 
showed that such models reproduce very well the ERO 
number counts in the range 18 ;$ K ^S 20, consistently 
with no appreciable evolution in number density up to 
z ~ 1.3. In Fig. we show some examples of the redshift 
selection functions of R — Ks > 5 ellipticals, as derived 
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Fig. 2. The selection functions of the ellipticals with R—Ks > 
5 for the passively evolving model described in the text are 
shown, for Zf — 2.5 (top) and Zf = 10 (bottom), and for 
various Ks limiting magnitudes, for the A-fiat cosmology. 



from our models, with various formation redshift and lim- 
iting Ks magnitude. 

In conclusion, we have modeled the selection function 
accounting only for passive evolution, since empirically the 
currently available best data including our own sample, in- 
dicate little (if any) density evolution and SED's and color 
trends consistent with pure passive evolution. Some uncer- 
tainty exists at z > 1.5 where negative density evolution 
could reduce the numbers of red objects in a magnitude 
limited sample, if merging has been important. However, 
the selection functions we obtain by setting Zf from 2.5 
to 30 bracket a large range of possible models and the dif- 
ference between these in terms of the high-z tail shown in 
Fig. and used in recovering ro, may be reasonably ex- 
pected to accommodate any modest density evolution like 
that claimed in the models of Kauffmann et al. (1999), 
Somerville et al. (2001). In fact, there are two features of 
the selection function that mostly influence the estimate of 
ro, i.e. the width Az and the effective redshift z e g (see Fig. 
g). Actually, the two limiting cases of Zf — 2.5 or Zf = 10 
shown in Fig. [2| produce only a modest ~10% variation in 
the ro estimate, basically because the two effects conspire 
to cancel each other (see Table |l|). 

4. Estimating the spatial correlation length with 
Limber's equation 

In D2000, clustering amplitudes for EROs with R—Ks > 5 
were estimated at several Ks limiting magnitudes (Fig. 0, 
see also Tab. 5 of D2000). By using the passive evolution 
dN/dz distribution appropriate for the different Ks limit, 
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Fig. 3. The data are the angular clustering measurements 
taken from D2000 (Tab. 5, filled circles). The empty square 
is the McCarthy et al. (2000) measurement with H converted 
to Ks using H — Ks — 1. The passive evolution models predic- 
tions are also shown, in the case of a A-flat cosmology. The ro 
value adopted for each model is the best fitting one, as shown in 
Tab. QJ. The general trend shown here is unchanged in different 
cosmologies. 

Table 1. Real space correlation lengths for EROs, derived 
through Limber's equation, assuming the selection func- 
tions expected for the ellipticals in the passive evolution 
case. The correlation lengths ro are expressed in comoving 
hr x Mpc, but a proper scaling to h values different from 
those used in the models would require a recalculation of 
the selection functions. 
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the predictions for the angular clustering amplitude have 
been derived at the same Ks limits by means of eq. 0, 
as a function of the ro values. The best estimate for ro 
was than obtained from a \ 2 minimization between all 
the predicted and observed angular clustering amplitude 

(Fig.|). 

Tab. [l] reports the inferred results, for different cosmo- 
logical models and redshift of formation. The Table shows 
that for any given cosmology the best fit values for ro are 
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Fig. 4. The variation of the effective redshift z e a (defined in 
eq. 0) for samples of passively evolving ellipticals, selected with 
R — Ks > 5, as a function of the Ks limiting magnitude, for 
the A-fiat cosmology. 



not a strong function of the unknown formation redshift. 
For each cosmology the worst agreement, as judged from 
the x 2 values (see also Fig. 0) is obtained with the small- 
est value of Zf . For each entry in the Table we estimate a 
typical error of the order of Aro ^S 1 h" 1 Mpc, obtained 
by propagating the measured AA values through eq. 0, 
for the three single most precise A measurements (the er- 
ror corresponding to a A\ 2 — 1 variation are significantly 
smaller than that). Given the small variations of r$ in 
Tab. pi as deduced with different formation redshifts, and 
their internal variance, we can conclude that, according to 
the Limber equation, EROs have a comoving correlation 
length of ro ~ 13.8 ± 1.5 h^ 1 Mpc in the A-fiat case, or 
ro ~ 11.5 ± 1 in the open or fi-flat case, applying to an 
effective redshift of 1 ,$ z e g ,6 1.2. 

Fig. H shows that the dependence of A on the Ks lim- 
iting magnitude is consistently reproduced by the dN/dz 
variations, so that the effective correlation length is not 
very dependent on the Ks magnitude within the observed 
range. This reflects the expected small variation of z e & 
over our samples (from z e g ~ 1 at Ks = 18 to z c g ^1.2 
at Ks — 19.2, see Fig. 0). We stress that the optimal 
agreement between the predicted and observed trend of 
the ERO angular clustering versus limiting Ks magnitude 
is good evidence of the consistency of our modeling of the 
selection function based on the passive evolution of the 
stellar populations. 

We also plot in Fig. K3 the angular clustering measure- 
ment by McCarthy et al. (2000), converting their limit of 
H < 20.5 by using the typical color of EROs H — Ks ~ 1 
(Cimatti et al. 1999), and assuming an error on their mea- 



surement similar to our best ones, given their total area of 
1000 arcmin 2 . Their color selection criterion of / — H > 3 
is roughly consistent with our R — Ks > 5. The McCarthy 
et al. (2000) point is in good agreement, at least within 
~ 1<t, with our model's predictions for Zf > 4 and with 
the general trend of amplitude versus limiting magnitude 
inferred by the D2000 survey. 

5. Estimating the spatial correlation length from 
numerical simulation 

Can we trust Limber's equation results, given the known 
spikey structure of redshift distributions in pencil-beam 
surveys? The size of our field (22 x 32 arcmin) corresponds, 
in fact, to ~ 17x 25 h~ 1 Mpc at z — 1, while in the redshift 
direction we sample objects over a range of approximately 
1000 hr 1 Mpc (A-fiat cosmology). The comoving correla- 
tion lengths derived from the Limber equation analysis is 
therefore of the same order of the projected size of our sur- 
vey, suggesting that a very clumpy redshift distributions 
should be expected for our sample. As Limber's equation 
has formally no dependence on the extension of the field 
over which the angular correlations are measured, we ex- 
pect that it should apply in the limit in which the analyzed 
field of view is large enough so that the observed redshift 
distribution approaches the selection function. But it is 
not clear a priori if the Limber's equation should apply 
accurately for small fields of view, in which the sampling 
of redshifts along the line of sight will vary considerably, 
dominated by notable spikes (e.g. Broadhurst et al. 1990, 
Cohen et al. 1999, Yoshida et al. 2001). 

Since this point has not been investigated previously 
we have embarked upon detailed modeling with realistic 
simulations. To do this we created clustered 3D distribu- 
tions of objects with known input ro in a suitable comov- 
ing volume, we then applied the selection function dN/dz 
to these simulated data and projected them on the sky 
for directly measuring the two-point angular correlation 
function for comparison with the data. 

To build up the 3D clustered samples the Soneira & 
Peebles (1977, 1978) algorithm was used with a 15 level 
hierarchy, setting the first pass radius equal to 50 h~ x 
Mpc and the step distance of the algorithm to the proper 
value to obtain correlations with 7 = 1.8. We refer to the 
original papers for details and discussions about the algo- 
rithm. By directly measuring the 3D two-point functions 
we calibrate the algorithm's parameters in order to repro- 
duce the desired ro value. Such measurements were done 
by a simple generalizations to 3D, following the approach 
summarized in D2000. The precision we reach is better 
than ro/Aro ^ 50 over the range 6 ^S r$ ,$ 15. 

The simulations were aimed at reproducing the ob- 
servations for the EROs with Ks < 18.8 for which we 
estimate an angular correlation of A = 0.014 ± 0.002, as 
this uses the full 701 arcmin 2 area (the measurement at 
Ks = 19.2 has a better S/N but is limited to a smaller 
area of 447.5 arcmin 2 ). Given that the analysis based on 
Limber's equation show that different Zf and different cos- 
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Fig. 5. Examples of simulated realizations of our field of 22 x 32 arcmin populated with a population with fixed 3D clustering 
ro = 12 h^ 1 Mpc. For reference, the top-right panel show the sky distribution of our data. 



mologies yield very similar results, we restricted our sim- 
ulations to the case of dN/dz produced by the model with 
Zt = 4 and adopted a A-flat cosmology. The idea is to 
build up a test case to better understand the behavior of 
projected clustering since we do not expect this to depend 
significantly on the details of the selection function. 

We have produced 120 independent realizations of a 
field of view of 22 x 32 arcmin to match the data, for each 
of 33 values of ro ranging from about 6 h~ x Mpc to 15 h~ l 
Mpc, and then we have measured the angular clustering 
in the same way as for the data (basing on the Landy 



& Szalay 1993 estimator, see D2000). The simulations are 
populated in such a way to produce on average 280 objects 
for each field to match the data for the EROs with Ks < 
18.8. In Fig. Hand |we show some examples of simulated 
sky projections together with their corresponding redshift 
distributions, for a population with ro = 12 h~ l Mpc (that 
we show in section 5.3 to be the best fitting value for 



EROs). 

The main results of the simulations are shown if Fig. 
M. These simulations show two interesting findings. First, 
the presence of a large dispersion in the measured values 
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Fig. 6. Examples of redshift distribution recovered from the 
simulations. Each panel refers to the corresponding one in 
Fig. B. In the top-right panel the adopted selection function 
is shown. 



of the amplitude of the correlation function, for any given 
assumed tq value. This in turns implies, for any given mea- 
sured amplitude, a large range for the statistically accept- 
able values for tq, well in excess of what obtained from 
the Limber's equation (see previous Section). Secondly, 
for each rg value the mean clustering amplitude recovered 
from the simulations is systematically higher than what 
predicted by the Limber's equation. 

5.1. Intrinsic variance of the two-point correlation 
function 

The origin of both these effects (large dispersion and bias 
of the amplitude of the correlation function) can be esti- 
mated from simple considerations. The basic line of rea- 
soning is the well known fact that, because of clustering, 
the actual variance in the object number counts is larger 
than the poissonian variance, and can be written as (see 
D2000, eq. 8): 



n (1 + nAC) 



(8) 



where n is the mean expected number of objects and AC 
is the integral constraint (Groth & Peebles 1977): 



AC 



1 

rJ2 



w{0)dn 1 dn 2 



(9) 



The expression for the variance in eq. M is the same 
which one would obtain if all the observed objects would 
belong to clumps with: 
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Fig. 7. The data points show the mean angular clustering am- 
plitude recovered from the simulations, as a function of the in- 
put tq. The error bars correspond to the standard deviation of 
the distribution of the recovered amplitudes. For comparison, 
the predictions of the Limber's equation are shown (curved 
line) . The shadowed area brackets our observed amplitude and 
its uncertainty A = 0.014 ± 0.002. 



(i.e., from eq. 0, the number of clumps is the inverse of 
the average of w(0) over the observed field), and the num- 
ber of objects per clump were a stochastic variable with 
average n/N c \. From eq. [To] it is expected that a variance 
in the number of clumps N c \ should result in a variance 
in the clustering amplitude A. In the minimal hypothesis 
that the clumps are distributed at random in the sky (i.e. 
neglecting the clustering between clumps) then: 



a N cl 



dN ci 
dA 



<*A 



(11) 



vvci = (Acy 



(10) 



from which it follows that the relative dispersion on the 
amplitude of the correlation function caused by the sky 
fluctuation of N c i would be: 



(12) 



Thus, one should observe a real variation of the corre- 
lation amplitude A on the sky, even at fixed 3D clustering 
length. The observations of angular clustering for a popu- 
lation of fixed r should result in a distribution of values 
with a variance decreasing with the area over which the 
measurements are carried out (the factor C is a decreas- 
ing function of the area, see eq. 9 in D2000), and strongly 
increasing with the expected average angular amplitude 
(<JA oc A 3 / 2 ). 

In any generic angular clustering measurement such 
intrinsic variance, that depends on the survey geometry 
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Fig. 8. The filled squares show the standard deviation versus 
the mean angular amplitude recovered from the simulations. 
Empty squares are the average statistical uncertainties in the 
A measurements. The solid line corresponds to the predictions 
of eq. fly, while dashed line is obtained by adding in quadrature 
both sources of variance. 



and the clustering strength, has to be added to the statis- 
tical uncertainty in the measurement of A that is linked to 
the finite number of observed galaxies. In principle with 
very large areas (and/or weak clustering) only the latter 
has a measurable effect, but in the case of small fields, if 
the clustering itself is strong, the former may become the 
dominant source of uncertainty. 

Our analytical derivation is supported by and can ex- 
plain the results of the simulations that we have carried 
out. Fig. shows that only by adding both contributions, 
the variance measured in the simulations can be recov- 
ered rather well, with some underestimation ( ;$ 10%) for 
large A values. Probably such small underestimation is 
explained by other effects neglected here, such as for in- 
stance the variance in the mean redshift of the clumps in 
a given beam, as the angular clustering is increased for a 
lower mean redshift. Incidentally, we note from Fig. |8|tliat 
the statistical uncertainty in each single measurement of 
A is found to be rather constant as a function of A, thus 
depending only on the number of observed objects. At the 
same time such statistical uncertainty seems to be overes- 
timated by a factor of ~ 2 (we recall that it follows from 
assuming 2 sigma poisson errors for the correlations as 
suggested by bootstrap analysis, see the discussion about 
it in D2000), given that lim^^o a A is about half of the 
average statistical uncertainty in the measurements of A 
(see Fig. |). 

The basic results of this analysis is that the variance in 
the angular clustering can be much larger than the purely 
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Fig. 9. The same of Fig. |, but for a field of 7000 arcmin 2 



statistical variance, especially for small fields of view. A 
similar qualitative conclusion was empirically reached by 
Postman et al. (1998) by splitting their large photometric 
survey into 250 smaller subunits, finding a large scatter in 
the amplitudes recovered from the smaller fields; a point 
taken up by McCracken et al. (2000) in the analysis of 
a single very deep field of 50 arcmin 2 . For the first time 
here we quantify this effect and give a general analyti- 
cal prescription to predict its amplitude. In the literature 
this additional source of variance is usually not considered, 
while deprojection analyses have been carried out even for 
surveys covering tiny fields of view of only a few arcmin 2 
of sky (e.g. from the Hubble Deep Fields) which must 
severely underestimate the true variance if the standard 
Limber's based inversion procedure is applied. Moreover, 
as a large variance is indeed expected for the measure- 
ments of A, our findings may help to explain why so many 
apparently discrepant clustering measurements (for simi- 
lar observing conditions) have been found in the literature 
(see e.g. Fynbo, Freudling & Moller 2000, McCracken et 
al. 2000). 

5.2. A possible bias in the Limber's equation based 
inversion 

The other interesting finding of our simulations is the pos- 
sible presence of a small (about 15%) but significant bias, 
with respect to Limber's equation predictions on the re- 
lation between tq and the measured angular correlation 
amplitude (Fig. |^). 

Also this effect can be understood by thinking in terms 
of objects coming in clumps, with A and N c \ linked by eq. 

because of 



nO. In fact, since in general < -^- >^ 



<A r d>' 
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Fig. 10. The probability to produce an ERO angular ampli- 
tude A > 0.014 as a function of the correlation length ro. The 
dashed lines show the la range for ro. 



eq. M we should expect deviations of < A > from the 
Limber's equation predictions if N c \ is small or if the dis- 
tributions of N c \ is significantly skewed. The inverse of 
A (proportional to N c \) should instead be less affected, 
and in fact we find that the average of the inverse of A 
is in better agreement with Limber's equation predictions. 
This effect is anyway small compared to the huge variance 
and so for our sample it is not an important correction. 
We have examined the possibility that this bias is at least 
in part produced by our calculation procedures generat- 
ing some systematic effect, such for instance a bias in our 
estimator (Hewett 1982; Kerscher et al. 2000) or an over- 
correction for the integral constraint. To rule out this pos- 
sibility we built up uncorrelated (random) samples with 
the same average number of objects and measured their 
two point angular correlations. This results in an ampli- 
tude of A = 1.4 10~ 4 ±1.2 10~ 3 , thus showing that pos- 
sible systematic effects on the estimation of A are much 
smaller than the bias we find. To verify further this point 
we carried out simulations as in Sect. pL over a larger area 
of 7000 arcmin 2 area (i.e. 2 square degrees, 10 times larger 
than our ERO field) keeping the surface density of objects 
fixed at the observed level. This exercise confirms the pres- 
ence of some bias although, as expected, at a lower level, 
and still increasing with increasing r suggesting that it is 
a real effect (Fig. M). This larger simulation also allowed 
us to verify that the trend predicted by eq. O still holds 
correctly. 



5.3. Application of the simulations to EROs 

From Fig. M it can be seen that the ERO observed cluster- 
ing amplitude A = 0.014±0.002 corresponds to ro ~ 12±3 
h^ 1 Mpc. This is considerably different from the simple 
Limber's equation estimate of ro ~ 13.2 ± 0.8 which we 
infer underestimates the uncertainty in the clustering am- 
plitude by a factor of ~ 3. Thus, the cosmic variance is 
a substantial source of uncertainty in the ro estimate for 
EROs, reducing the weight of the possible uncertainties 
in our selection function modeling. Using the spread in 
the A measurements from the simulations shown in fig. 
[7| we may place a lower limit to the correlation length of 
ro it 8 h^ 1 Mpc (see Fig. |l(j) at the 95% confidence level. 
While this estimate is derived only from the measurement 
at Ks = 18.8, we expect it to be consistent with what we 
would have found from the analysis of the clustering at all 
Ks levels. This is because within each single survey the 
number (and the redshift) of the clumps is fixed in the real 
space, so that we expect that the trend of the clustering 
amplitude as a function of the apparent magnitude should 
basically reflect only the change of the redshift distribu- 
tion (see Fig. ||). 

Our simulations have been carried out for convenience 
for the A-flat universe. From the results of Sec. we can 
say that in the open or fi— flat cases the ro estimate is 
lower by a small amount, so that at the 95% confidence 
level it becomes ro ^ 7 hr 1 Mpc. 

Fig. H shows that the angular clustering measurement 
of McCarthy et al. (2000) is consistent with our measure- 
ments, even if slightly lower than our model predictions. 
Given the large variance expected, the discrepancy is not 
significant. Nevertheless McCarthy et al. (2000) estimate 
from their data ro ~ 8 hr x comoving Mpc, significantly 
lower than our preferred ro value. The main reason for 
this lower ro value is in their adoption of a relatively nar- 
row Gaussian form of dN/dz, centered at z — 1.2 and 
with <j z = 0.15, motivated by their photometric redshift 
estimates. Such a strong confinement of EROs into a nar- 
row redshift range is not in agreement with our modeling 
of the selection function shown in Fig. |[ Even if Fig. g 
demonstrates that occasionally in small fields the observed 
dN/dz could be spuriously narrow (because of the cluster- 
ing) , the inversion process requires the use of the actual se- 
lection function, which we expect to be much broader than 
the observed, clumpy dN/dz of a given field. Assuming 
that the McCarthy et al. threshold of / — H > 3 is consis- 
tent with R — Ks > 5, and using our estimate of the ERO 
selection function at Ks < 19.5, the McCarthy et al. an- 
gular measurement could be inverted to ro = 10.8 ± 2.2, 
where the uncertainty is derived from our own one by 
keeping into account the scaling with the area and clus- 
tering amplitude, consistently with eq. [l^. This estimate 
apply to an effective redshift z e s ~ 1.2 (Fig. Eh. 
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6.1. Comparison to the clustering of bright local 
ellipticals 

We now compare the large correlation length estimated for 
the z Z, 1 ellipticals with that of their local counterparts, 
in order to constrain the cosmic evolution of the clustering 
of massive early-type galaxies. 

The correlation length of a population of galaxies is 
known to depend on the absolute luminosity selection 
threshold, and can also depend on the scales over which 
the clustering is measured. Such quantities must be prop- 
erly estimated in order to compare the clustering of EROs 
to that of local ellipticals. For a passively evolving ellip- 
tical, the apparent magnitude of Ks — 19.2 corresponds 
to L ~ 0.6L* and L ~ 1.31* at the redshift of 1 and 1.5, 
respectively, while for Ks = 18 we have L ~ 1.6L* at 
z — 1 and L ~ AL* at z = 1.5 (accounting for the passive 
evolution of L*). Therefore our sample consists of galaxies 
with typical luminosity L £, L*. The largest effective sep- 
aration probed by our angular clustering measurements is 
8 ~ 15', corresponding to about 12 hT 1 Mpc at z ~ 1 
(A-flat universe). 

As the clustering amplitude is expressed with j-q, the 
measurements of r$ obtained with a 7 different from the 
value adopted here must be rescaled to that value in order 
to produce a meaningful comparison. Therefore, all the r 
quoted below were transformed with 7 = 1.8. 

Our results can be compared with those obtained lo- 
cally, for two different samples, by Guzzo et al. (1997, the 
Perseus Pisces redshift survey) and Willmer et al. (1998, 
the SSRS2 redshift survey). Both estimate the clustering 
of bright early type galaxies with Mb < —19.5 + 5logh, 
corresponding to L £ L*. Guzzo et al. measure scales up 
to about 10 hr 1 Mpc and find r = 11.3 ± 1.3 hr 1 Mpc, 
while Willmer et al. measure up to about 20 hT 1 Mpc 
and find tq = 7.6 ± 1.2. The two measurements are only 
consistent with each other at the 2<j level, but it should 
be noted that the Perseus Pisces redshift survey has a 
higher abundance of local clusters. It may be implied by 
Fig. 10 of Willmer et al. 1998 that they would obtain a 
larger amplitude if limited to smaller separations. For fur- 
ther constraints we note that local radio galaxies, known 
to be hosted by bright ellipticals, have r$ — 11 ± 1.2 hr 1 
Mpc (Peacock & Nicholson 1991). Therefore a correlation 
length in the range r§ ~ 9-11 hr 1 Mpc can be regarded 
as a reasonable estimate for the clustering amplitude of 
local L ^ L* ellipticals. 

Our estimate of ro = 12 ± 3 Yr 1 Mpc then implies 
that the clustering amplitude of bright ellipticals does 
not significantly decrease from to z ~ 1 — 1.5 to the 
present. The stable clustering scenario (i.e. e = 0, if 
r (z) = r (z = 0)(1 + z )(7-3- e )/7) that is k nown to g t 

many observations of clustering evolution to z ~ 1 (e.g. 
Peebles 1980, Le Fevre et al. 1996, Carlberg et al. 1997), 
would predict ro ~ 6 hr 1 Mpc at our inferred effective red- 
shift of Zed ~ 1.1 and hence it is not in good agreement 
with our measurement of the ERO clustering. A negative 
value for e is supported by our analysis. 



6.2. Comparison to theoretical predictions 

The evolution of the correlation function is popularly char- 
acterized as: 



az,r)=D 2 (z)b 2 (z)Ua SS (0,r) 



(13) 



where the purely linear growth D{z) of density pertur- 
bations, Ap(z)/p(z) is separated from the bias evolution 
b(z). In the linear case the bias evolution can be expressed 
with the Tegmark & Peebles (1998) prescriptions and well 
known expressions for the linear growth factor can be as- 
sumed (Peebles 1980, Treyer & Lahav 1996). In our case 
the linear assumptions are not likely to be correct, as for 
EROs we are mapping angular comoving scales similar to 
the inferred ro, thus sampling a region with £ ~ 1, so 
that this prescription can be considered as a rough bench- 
mark picture in the absence of the complexities affecting 
small scale growth (Kauffmann et al. 1999, Somerville et 
al. 2001). 

The simplest and most clear model of bias evolution 
of ellipticals is provided by the so called galaxy conserva- 
tion scenario (Fry et al. 1996, Tegmark & Peebles 1998, 
Moscardini et al. 1998, Magliocchetti et al. 1999, Lacy 
2000), that holds if the galaxy population is conserved over 
the cosmic time (i.e. no new elliptical forms and no one dis- 
appears). This scenario implies the assumption that all the 
ellipticals were formed at high redshift and simply follow 
the growth of perturbations without the additional non- 
linear effects such as virial collapse and merging expected 
at small scales (below ro of the mass auto-correlation func- 
tion). It is therefore relevant as a limiting case for compar- 
ison with our observations. Here the positive evolution of 
the bias which increases with redshift is more than com- 
pensated by the decline of linear growth, i.e. D{z) wins 
and a net decline of the clustering amplitude with increas- 
ing redshift is expected in the linear regime beyond ro of 
the mass-autocorrelation function. 

In Fig. [ll] (A-flat case) we see that normalizing the 
linear predictions to the ro and bias of local ellipticals 
(such bias was computed assuming ro = 10 h~ x Mpc and 
the normalization a™ ass — 0.9, from Eke Cole & Frenk 
1996), the predicted trend is slightly decreasing in comov- 
ing units, reaching values around ro ~ 7-8 h~ x Mpc at 
z = 1.1. Given the large uncertainties, this scenario can- 
not be rejected with great confidence, but it is disfavored 
by our findings. 

If the galaxy conservation model predictions, both for 
the bias and for ro, are extended to higher z, they do 
intercept the measurements for the LBGs clustering by 
Giavalisco et al. (1998) and Adelberger et al. (1998) (see 
Fig. O for the A-flat case). This kind of argument had 
led to previous suggestions that LBG's could evolve into 
present-day bright ellipticals (e.g. Adelberger 2000). An 
important confirmation of this picture would be to find 
that distant (i.e. z ~ 1) ellipticals have intermediate clus- 
tering strength between the local ones and the LBGs. Our 
findings, taken at face value, disfavor such an interpreta- 
tion as they seem to suggest that the clustering of chip- 



E. Daddi et al.: Spatial Clustering of z ~ 1 Early- Type Galaxies 



11 




Fig. 11. The ro measurements for local ellipticals (filled 
squares), z ~ 1 ellipticals (EROs; filled circle from D2000; 
the empty square is derived from our analysis applied to the 
McCarthy et al. angular clustering measurement) and LBGs 
(asterisks) are shown for the A-flat universe (see the text for 
details). The dot-dashed line shows the Kauffmann et al. (1999) 
ACDM predictions for the clustering of ellipticals. The case of 
stable clustering (e = 0) is indicated by the dotted line. The 
other three curves show the predictions for the galaxy conser- 
vation scenario, differing in the assumed degree of correlation 
between galaxy and dark matter, assumed from top to bottom 
to be 0.9-0.95-1 at z = (see Tegmark and Peebles 1998, for 
more details on this parameter). All the models were normal- 
ized to r (0) = 10 h" 1 Mpc. 



ticals is not decreasing enough from z ~ to z ~ 1 to 
reproduce the clustering at z — 3, at least for the LBGs 
with absolute luminosity similar to those in the Giavalisco 
et al. (1998) and Adelberger et al. (1998) samples. 

Alternative more complex models for the clustering 
evolution have been constructed to deal with nonlinear 
effects and may be termed the galaxy merging scenario for 
bias evolution (e.g. Mo & White 1996, Moscardini et al. 
1998), in which the bias grows with z with some law of the 
kind A6 ex (1 + z) 1 - 8 , which is stronger than the growth of 
perturbations, with the net effect that ro increase with z. 
Such models are in better agreement with the ERO clus- 
tering, especially if the large favored ro ~ 12 h^ 1 Mpc will 
be confirmed. Similar predictions for the clustering evolu- 
tion of the ellipticals in the ACDM semianalytic model of 
Kauffmann et al. (1999) and Somerville et al. (2001) are 
also consistent with the trend inferred here (see Fig. Ill]) 
and include further refinements such as luminosity evolu- 
tion. 

A difficulty one might expect of strong bias evolution 
models is that they may contradict the observational ev- 
idence of the lack of evolution of the comoving density of 



ellipticals, which appears not to decrease significantly to 
z = 1 and beyond (see Daddi Cimatti & Renzini 2000). 
However, depending on the details of the semi-analytical 
approach, it is possible to accommodate substantial bias 
evolution of the dark matter without perturbing either the 
apparent space density or clustering amplitude of ellipti- 
cal galaxies for the ACDM model (Kauffmann et al. 1999; 
Somerville et al. 2001, see also Bullock et al. 2001) and the 
main difference between these models and the conserva- 
tion scenario would seem to be the inclusion of a suitable 
prescription for merging. As discussed in Daddi Cimatti & 
Renzini (2000), the small amount of density evolution re- 
quired by the ACDM models could be consistent with the 
ERO number counts once it is required that merging does 
not produce significant star formation (that would make 
the objects bluer than our color cuts), i.e. red merging is 
required. This indeed is claimed to have been observed in 
clusters (van Dokkum et al. 2001). 

7. Summary and conclusions 

The main results presented in this paper are: 

• The real space correlation length of EROs (with R — 
Ks > 5) is much larger that the correlation of z ~ 1 star- 
forming galaxies for any reasonably large selection func- 
tion, strengthening the previous suggestions that most 
EROs at Ks ~ 19 are early-type galaxies. 

• Assuming that EROs are predominantly early-type 
galaxies and hence that their selection function is rea- 
sonably described by passive evolution, then the spatial 
correlation length we obtain is rather large, not less than 
7-8 hr 1 Mpc, with the most probable estimate of 12 ± 3 
h^ 1 Mpc, applying to an effective redshift of z ~ 1.2. 

• At face value this implies no significant evolution of clus- 
tering of this population relative to the present day when 
measured in comoving coordinates, and a strong bias in- 
crease from z = to 1. 

• Realistic simulations were used to constrain ro from the 
observed angular clustering of EROs. Two main results 
follow from the simulations that can be of general inter- 
est for the analysis of the angular clustering: (1) the am- 
plitude A of the angular two-point correlation function 
fluctuates on the sky according to a a/ A = VAC, and 
(2) a possible systematic overestimate of r could follow 
from the inversion of the angular clustering measurements 
based on Limber equation. Both effects are strongly en- 
hanced in the limits of strong angular clustering and/or 
small fields of view. 

Taken at face value, our result on the ERO correlation 
length challenges the simple conservation models of clus- 
tering growth for massive haloes, but it may be reconciled 
with more complex schemes for the bias which incorpo- 
rate more parameters to describe non-linear effects such 
as merging, so that although evolution of the underlying 
mass function strongly declines with redshift, the observer 
will, it is claimed, find the opposite of the expected behav- 
ior, namely an increase in the observed correlation length 
of early-type galaxies with redshift and no corresponding 
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strong decline in their number density with increasing red- 
shift (Kauffmann et al. 1999, Somerville et al. 2001). We 
rule out a high degree of density evolution of early type 
galaxies to z ~ 1. This is inconsistent with the ERO clus- 
tering because strong density evolution would significantly 
narrow the width of the selection function by removing the 
high-z tail resulting in a reduction of the inferred estimate 
of ro. This is counter to the strongly increasing correla- 
tion amplitude that would be expected with increasing 
redshift for such a highly biased model. A modest reduc- 
tion in density at z > 1 can be accommodated given the 
present uncertainties in lookback time and star formation 
which fold into the construction of the selection function. 
Indeed some change in density through merging is sug- 
gested by our results when we combine the constraints on 
both density and correlation length evolution. 

Before discarding the galaxy conservation scenario for 
the clustering evolution of early type galaxies some possi- 
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